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Abstract. We have studied the Z2 topological number defined by the Berry phase for the 
gapped frustrated systems including the orthogonal dimer model which has a direct product 
state of local quantum clusters as the exact ground state. The topological number can 
clarify what kind of the local quantum clusters is formed to lift the macroscopic degeneracy 
due to frustration, even when the exact ground state is unknown. As a demonstration, the 
dimer-singlet and plaquette-singlet phase are identified by two kinds of Z2 topological numbers 
in the Shastry-Sutherland model and its generalization realized experimentally as SrCu2 (603)2 
and CaV4 09. 



1. Introduction 

Macroscopic degeneracy that survives at low temperature is a key feature of frustrated materials. 
Approaching the absolute zero temperature, the degeneracy will be lifted to meet the third 
law of thermodynamics. How is the degeneracy resolved? Its answer has a rich diversity, 
because frustrated systems have many exotic phases that are sensitive to small change. A 
standard scenario is the symmetry breaking mechanism that tells us an order parameter to 
characterize how the degeneracy is resolved. However, recent development in condensed matter 
physics reveals that there are topological phases where usual order parameters do not play 
any fundamental roles. Instead of usual order parameters such as a magnetization, topological 
invariants are used to characterize the topological phases and they are linked to the presence 
or absence of gapless edge modes. A historical example is the quantum Hall system where the 
topological invariant is the Chern number that is an integer defined by the Berry connection 
theoretically and observed as the Hall conductance experimentally. In addition to the Chern 
number, Z2 topological values, i.e., the integers modulo two, have been focused in recent studies 
on quantum spin Hall systems[l, 2]. In such topological phases, many of them are gapped. 

In some frustrated materials, gap formation is another possible scenario for lifting the 
macroscopic degeneracy as in the two-dimensional gapped spin system, SrCu2(B03)2[3, 4]. A 
discovery of SrCu2(B03)2 shed a light on theoretical construction of the model Hamiltonian with 
an exact ground state of local spin-singlets, i.e., the Shastry-Sutherland model[5, 6]. Inspired by 
Majumdar's one-dimensional model[7], Shastry and Sutherland constructed the two-dimensional 



Figure 1. Gap AE' and Berry 
phases of the Shastry-Sutherland 
model for system size = 16 with 
the adiabatic parameter a = J'/ J. 
J (and J') is the nearest neighbor 
(next nearest neighbor) interaction. 
7d and 7p are Berry phases for the 
DS and PS, respectively. The detail 
is discussed in §3. 

geometrically frustrated spin system in which the exact ground state is the dimer-singlet(DS) 
state which is a direct product state of local spin-singlets [5]. Unless the spin gap closes, the 
DS state is still the exact ground state against the particular kind of perturbation due to the 
orthogonality of the dimers. Such construction is common in orthogonal dimer models including 
three-dimensional generalization [8] . 

The Shastry-Sutherland model also has an gapped phase, the plaquette-singlet(PS) phase[9- 
11], in addition to the DS phase. Unlike the DS phase, it is difficult to obtain its exact ground 
state in the PS phase. In fact, Weihong eta/. [8] have raised a possibility that the PS phase is 
unstable. To clarify it, one can consider an adiabatic modification of the Hamiltonian connected 
to the decoupled PS Hamiltonian which has a direct product state of PSs as the exact ground 
state. Actually, the ground state in the PS phase is adiabatically connected to that of the 
decoupled PS Hamiltonian [9]. The generalized theoretical model discussed in the adiabatic 
connection includes the 1/5-depleted square lattice that is realized as CaV409[12, 13]. In this 
sense, the direct product state of local quantum objects, such as DSs and PSs, and the adiabatic 
modification toward decoupled Hamiltonians are not just theoretical toys. 

Recently, a method to obtain a decoupled Hamiltonian is proposed, that is, the Z2 topological 
number defined by the Berry phase[14-18]. Unless the spin gap closes, the Z2 topological number 
remains invariant, i.e., topologically protected against perturbation. These situations are quite 
similar to the case of the Chern number. The difference is that the Chern number is integer 
while the Z2 topological number is quantized to the following two values: a trivial value or 
nontrivial value vr (mod 2ti) due to the time reversal symmetry of the spin system. However, 
since the Berry phase is generally defined through the local perturbation, we can consider many 
kinds of the Z2 topological numbers defined by the Berry phase. The Z2 topological numbers 
have successfully identified the spin-singlets, plaquette-singlets, and Haldane state in the spin 
ladders. It can be applied to not only the spin system but also the electron system including 
the BCS system[18]. In addition, the itinerant singlets and the Kondo singlets are identified. 
Recently, we have succeeded in generalizing the Z2 Berry phase into Zq Berry phase with any 
integer (5[19]. 

In this paper, we apply the Z2 topological number to the Shastry-Sutherland model. Our 
numerical result in Fig. 1 shows that the DS and PS phase are clearly identified by two kinds 
of Z2 topological numbers, 7d and 7p, respectively. Numerical detail will be explained in §3, 
but here we note that the ground state in the PS phase is not a simple direct product state, 
while there is the simple exact ground state in the DS phase. Generally speaking, although the 
direct product state is not an exact ground state when small long range Heisenberg exchange 
interactions are considered as a perturbation for real material of SrCu2(B03)2, it will be widely 
accepted that the ground state in the perturbation is similar to the direct product state. The 
similarity is naively based on the adiabatic connection of two ground states. By using the 
Z2 topological classification, it is possible to discuss clear equivalence of ground states via the 




adiabatic connection, as we will explain in §2. This approach is demonstrated for the generalized 
theoretical model including SrCu2(B03)2 and CaV409 in §3. 

2. Definition of the Z2 topological number for orthogonal dimer models 

To describe our approach, let us start with local quantum clusters. A local quantum state IV'i) 
is supposed to be the singlet ground state of a quantum local Hamiltonian hi of an ith cluster, 
such as a spin-singict, plaqucttc-singiet, local Kondo-singlet and dimcrizcd electrons state. The 
decoupled Hamiltonian is given by Hq = J^i hi with [hi, hj] = 0. Due to small system-size of hi, 
the energy gap A generally opens above the exact ground state |*o) = Hi IV'i)- Although the 
local singlet-cluster \ipi) is a pure quantum object and does not have a classical analogue, the 
direct product state |5'o) is a classical state in the sense that there is no quantum entanglement 
between the clusters. The decoupled Hamiltonian is a start point to consider both construction 
of the orthogonal dimer model and definition of the Z2 topological number. 

To construct an orthogonal dimer model generically, we connect the decoupled Hamiltonian 
Hq with an interaction Hi. Introducing an adiabatic parameter a, we define = (1 — a)Ho + 
aHi or Ha = Hq + aHi. To keep the direct product state l^'o) as the exact ground state for 
a > 0, a condition Hi\^q) = is required. Since we supposed is a singlct-state. Hi can be 
constructed by a total spin operator for ith cluster T.^ = J^i.Sfi {P = x,y,z) with T^\%lji) = 0, 
where il indicates the /th site in the ith cluster. The condition guarantees that the ground state 
of Ha is the direct product state around a = to some extent as far as the gap is non-zero. 

In the orthogonal dimer model, we set the local spin 1/2 cluster hi = JSn ■ Sj2, with the 
local singlet ground state = t)ii| 4-)j2 — I 4-)ii| t)j2)- Using the total spin operator 

Tf = + S^2, we can construct Hi = Ju'v^i'V • (Sa + Sj2)- The coefficients J|^,;, are 

determined by the lattice geometry. In this construction, the orthogonality of the dimer bonds is 
always existing, that is, the interaction between il and i'l' is identical to that between i2 and i'l' 
always. Hi can be generalized to a scalar product type Hi = J2 '-^ijpj'p' ^ ^j'P') ' + 
and more complicated types. In addition, the local quantum objects can be generalized to 
a plaquette singlet when we consider hi = J(Sji • Sj2 + 8^2 • 8^3 -|- 8^3 • 8^4 -|- 8^4 • 8^1) with 

= Y^i=iS^i = Tf\ In the same manner, one can consider larger cluster singlets and 
generalizations to local quantum objects of bosonic or fermionic systems. Note that the Hermite 
condition Hi = hI imposes some limitation on Hi. In the above construction, the gap A of Hq, 
which comes from a finite-size gap of the local quantum cluster, plays an important role for the 
stability of the exact ground state against the adiabatic parameter a. 

If the local ground state \ipi) is not a singlet but a multiplct, we must deal with the 
macroscopic degeneracy. The gap will open above the macroscopically degenerate states and is 
identified by the Z2 topological number defined by non-Abelian Berry phases. Of-course, the 
perturbation Hi may lift the degeneracy and we need further discussion if we focus on the states 
below newly opened gap. In the following, we limit ourselves to a non-degenerate ground state 
and Abelian Berry phases in the definition of the Z2 topological number. 

The Berry's geometrical phases [20] are used to calculate the macroscopic polarization in the 
solid with surfaces. We emphasize that an edge or surface property such as the macroscopic 
polarization can be obtained via the Berry phase which is calculated in the periodic bulk-system 
without surfaces. This is bulk-edge correspondence, which is originally proposed as the one- 
to-one correspondence between the Chern number defined by the Berry connections and the 
number of the edge states in the quantum Hall systems[21]. Such local degrees of freedom 
on the boundaries are also characteristic in the topological insulators, and quantum spin Hall 
systems [1, 2]. The bulk-edge correspondence is important for the Z2 topological number used in 
this paper. For example, the non-trivial Z2 topological number guarantees the existence of the 
Kennedy triplet for an open chain[17]. 



To illustrate the definition of the topological number for B.^-, let us consider the decoupled 
Hamiltonian first. In our previous studies, it have been clarified that the gapped local 
quantum cluster, such as the spin-singlet, plaquette-singlet, and Kondo singlet, gives the non- 
trivial TT Berry phase. The key ingredient to obtain the vr Berry phase for the local cluster 
hi is the unitary gauge transformation. Let us explain it with detailed definition of the Berry 
phasc[21]. The Abclian Berry phase is defined by a local SU(2) spin twist on a specified site. 
When we choose a local cluster hi and a Ith site in hj, the local spin twist is introduced by 
the substitution hi hi{(f)) = U^(f))hiU{(f)) with = exp [i(l){S - Sfi)] for a spin 5 = 1/2, 
where II indicates the chosen site. With using hi {(()), the decoupled Hamiltonian has one- 
parameter dependence as Ho{(p) = hi{(f)) + J2i'^i ^i'- With the adiabatic parameter a we define 
Hai(p) = (1 - a)Ho{4>) + aHi or = Ho{4>) + aH^. 

For an one-parameter dependent Hamiltonian Ha{4>), the Berry phase 7 is defined as 
1 ~ lo^ cl<A(^Q(0)l"^|'J'a(0))) where \'^a{(l>)) is the ground state of Ha{(t)). The numerical 
integration is given by discrctizing the parameter space of (j) into finite number of points[22], 
which is enough up to 60. The Berry phase here is quantized into 0, or tt if the ground state 
is invariant under an anti-unitary operator. The anti-unitary operator is the time-reversal 
operator for the spin system. For the decoupled Hamiltonian Hq{(/)), the ground state is a 
direct product state and has one-parameter dependence through the unitary transformation 
\^o{4>)) = U{(p)\'i'o{0)), which gives the tt Berry phase, because 7 = 27r(5-(^'(0)|S'f;|'I'(0))) = tt 
in the spin S* = 1/2 case. This result does not depend on the detail of the local Hamiltonian 
hi. When the Hamiltonian is Ha{(j)) = (1 — Ci)Ho{(f)) + aHi, there should be the change of 
the Z2 topological number. At a = 1, since the Hamiltonian i?Q=i(0) = Hi does not depend 
on 4>, the Berry phase 7 cannot become the non-trivial value tt, that is, 7 is zero if the gap 
opens or undefined if the gap closes. Then, there is inevitably the change of 7 in the adiabatic 
modification via a, which signals whether the ground state is equivalent to the direct product 
states or not. It should be noted that since Ha{4>) / {(t))HaU{(t>) due to i^i, the Hamiltonian 
Ha{4') is able to have a local flux, which cannot be gauged out. Then, the problem is non-trivial 
for generic Hi. 

Typically, there are two types of the transitions for the Z2 Berry phase in the adiabatic 
modification a. To discuss it, let us show an example of a local cluster, the N = A spin 
S = 1/2 system with the decoupled Hamiltonian HQ{(f} = U'^ {(l))hiU {(j)) + /12 of the dimer 
singlets hi = Sji • Si2 and the local spin twist U{(j)) = exp S* — 5"^^)]. Here we consider 
two kinds of Hamiltonians Ha{4') = (1 — a)i^o(0) + oiHi: (a) the non-orthogonal interaction 
Hi = S12 • S21 + Sii • S22 and (b) the orthogonal interaction Hi = (Sn -|- S12) • (S21 + 822)- Due 
to the lattice symmetry, the Berry phase does not depend on the site where U{4>) is introduced. 
In the case (b), the exact ground state of Ha{0) is a direct product state of the dimer singlets 
for a < ac and there is a level cross at a = ac = 1/2. On the contrary, in the case (a), the direct 
product state is not an exact ground state for a > and there is finite gap without any level 
cross in the range of a G [0,1] for Ha{0). We expect that in the latter adiabatic modification 
two DSs delocalize asymmetrically except for the symmetric point a = ac- The Z2 Berry phase 
detects the asymmetry of the DSs. In Fig. 2, the result of the Z2 Berry phase is shown in both 
cases (a) and (b). In the case (a), the Z2 Berry phase clearly identified the two phases connected 
to the two decoupled models, a = and 1, respectively. At the phase boundary a = ac, the 
Z2 Berry phase becomes undefined because the level cross occurs at = tt. In other words, 
the diabolic point of the Dirac cone generates the tt Berry phase. Here, we emphasize that the 
Hamiltonian Ha{0) at cp = has no gap closing in the adiabatic modification via a and there 
is not a transition but a cross-over [18]. In the case (b) as shown in Fig. 2(b), the level cross at 
a = ac without twist <p induces a finite range of the undefined region of the Z2 Berry phase. 
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Figure 2. Two typical 
energy diagrams as a function 
of (j) and a of Ha{(j)) = 
(1 - a)Ho{<p) + aHi with (a) 
the non-orthogonal interaction 
Hi = Si2 • S21 + Sii • S22 and 
(b) the orthogonal interaction 
Hi = (Sii + S12) • (S21 + 822)- 



3. Results on the Shastry-Sutherland model 

Here we discuss adiabatic modifications of the Shastry-Sutherland model as shown in Fig. 1. 
The Shastry-Sutherland model with J and J' = aJ can be defined as the case Ji = J and 
J2 = J3 = J' of the generalized model as shown in Fig. 3. As mentioned in previous sections, 
the definition of the Berry phase depends on the local spin twist and on the choice of the 
decoupled Hamiltonian for a given Hamiltonian. For this Shastry-Sutherland model, 7^ (and 
7p) is defined by the local spin twist cp^ (cpp) shown in Fig. 3. The decoupled DS (and PS) 
Hamiltonian only with Ji (J2) corresponds to 7^ (7p). Since the Shastry-Sutherland model has 
the exact ground state of the DSs, the level cross at a = ac is expected as in Fig. 2(b). Actually, 
in our calculation with the small system size A'^ = 16, there is the first-order quantum phase 
transition, that is, a gap closing A = occurs at ac = 0.67 as shown in Fig. 1. Moreover, 7^ 
shows a clear transition at a = ac in Fig. 1. This is quite natural because we know the exact DS 
ground state for a = J' /J < ac- The question arising here is how the ground state for a > ac 
is identified. To answer this question we introduce another Z2 topological number, 7p, which 
corresponds to the decoupled PS Hamiltonian having the direct product state of local PSs as 
the exact ground state. Figure 1 shows that the two Z2 topological numbers can identify the 
DS phase for a < Oc and the PS phase for a > ac- However, since there is the first-order phase 
transition at ac, the Z2 Berry phase will show instability just around the level cross point as 
shown in Fig. 2(b). In fact, 7p shows an undefined region around a = Oc as shown in Fig. 1. 
Although the system size is too small and the transition point ac is not comparable to previous 
studies, the Z2 Berry phase clarifies a PS-type character of the ground state for a > ac- 

In addition, we consider J2/J1 = 1-82 and J3/J1 = 0.95 in the parameter space Ji,J2, and 
J3 of the generalized model in Fig. 3, which corresponds to CaV409 and is in the PS phase[10]. 
To clarify the property of the ground state by the Z2 topological number, we calculated 7p and 
we obtained 7p = vr for A^ = 16 sites. To illustrate the result, we show the gap AE as a function 
of twist (p = (j)p and the adiabatic parameter a from the decoupled PS Hamiltonian at a = 
to CaV409 at a = 1. The decoupled PS Hamiltonian has 7p = vr by definition. Figure 4 shows 
that there is no gap closing in the adiabatic modification not only at i;^ = but in whole region 
of (j) € [0, 27r]. To change the Z2 topological number, the gap closing is required as discussed in 
§2. In this sense, 7p = vr means that the state is adiabatically connected to the direct product 
state in the two dimensional parameter space including the gauge twist angle cp £ [0, 2tt]. 

4. Summary 

In summary, we have studied the Z2 topological number of the orthogonal dimer model as 
the gapped frustrated spin system. The construction of the orthogonal dimer model with the 
exact ground state and the Z2 topological number 7d are based on the same decoupled DS 
Hamiltonian. As shown in Fig. 1, the DS phase has the non-trivial Z2 topological number, 
7d = TT, corresponding the direct product state of DSs. The orthogonal dimer model shows the 
first-order transition with a level crossing, as shown in Fig. 2(b), in the finite-size calculation. 
Although the exact ground state after the transition is unknown, the Z2 topological number 
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Figure 4. The gap AE for the Z2 topological number 7p. The 
twist (p = (pp is defined for the decoupled PS Hamiltonian with 
the bonds J2 = 1.82. The adiabatic parameter a is introduced as 
Ji = a, J3 = 0.95a. 



Figure 3. Hamiltonian 
with Ji , J2 , and J2 . 



which corresponds to the decoupled PS Hamiltonian shows jp = ir for a > Oc as shown in Fig. 1 
and for the specified parameter corresponding to CaV409 as clarified in Fig. 4. Even though our 
analysis using the exact diagonalization is limited for the small system-size, this result agrees 
with existence of the PS phase discussed in the previous studies[9, 10]. 
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